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(S) Abstract 

! , Present data prefer a large but non-maximal 2 — 3 mixing in the lepton 

sector. We argue that this value, in connection with sin #13 ~ 0.15, 
Qh is the generic outcome of minimal flavour structures. We present a 

O | few different incarnations of this statement, in terms of lepton mass 

D matrices depending on a small number of parameters, that can be 

justified by discrete flavour symmetries. We also propose a general 
procedure to study the correlation between 823, the absolute scale and 
t-H ordering of the neutrino masses, and the leptonic CP-violating phases. 

> 

1 Introduction 

CN 
in 

Flavour models aim to explain the observed pattern of fermion masses and 
mixing, in terms of some dynamical mechanism for flavour symmetry break- 
er ing. In the quark sector the regular hierarchy of masses and mixing angles 
definitely points to an underlying dynamics, that unfortunately may be dif- 
ficult to probe directly. In the lepton sector, in the last 15 years neutrino 
oscillation experiments led to a tremendous improvement in the knowledge 
of the neutrino masses and mixing. Still a few important observables in the 
neutrino mass matrix are presently unmeasured, and one expects substantial 
improvements in the next decade. 

A reasonable (but not exclusive) figure of merit for flavour models is the 
ability to predict a set of observables in terms of a smaller set of parameters. 
Speculations on "symmetric" values of the lepton observables (maximal or 
zero mixing angles, vanishing or degenerate mass eigenvalues, . . . ) have been 
pursued, and recent indications against these extreme values may seem to 
disfavour most theoretical proposals. However, one should keep in mind that 
minimal flavour structures do not imply extreme values of the observables, 
as we will show in a few explicit examples. 
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Oscillation experiments are sensitive, in particular, to the three lepton 
mixing angles. Before 2012, only $12 and 623 were known to be non-zero 
and large. One year later, the situation changed dramatically: the 1 — 3 
mixing angle #13 has been precisely measured by reactor experiments [I] to 
be of the order of the Cabibbo angle and, given the precise determination 
of #12 by solar experiments, the largest uncertainty pertains to 623, that 
could significantly deviate from the maximal value. Indeed, global fits [2j 
|3| favour non-maximal 2 — 3 mixing at the 2a level (an effect driven by 
accelerator neutrino data) , with a slight preference for 823 < 7r/4 in the case 
of normal ordering of the mass spectrum (as a consequence of atmospheric 
neutrino data). Before the measurement of a non-zero #13, a considerable 
theoretical effort was invested in explaining a tri-bi-maximal pattern for the 
lepton mixing angles (see Ref. [I] for a review). In the latter approach, one 
needs now to correct significantly the leading-order predictions, by taking 
into account a number of sub-leading effects. After the measurement of a 
relatively large #13, several flavour models have been proposed, that predict 
a non-zero value of #13 at leading order (5J [6] . 

In this paper we would like to reconsider the lepton flavour structure in 
the light of the up-to-date knowledge of the mixing angles. Our focus will be 
on the non-maximal value of #23 and on the search for the minimal flavour 
patterns that lead to such value. Rather than elaborating on general aspects 
of model-building with flavour symmetries, we will modestly propose a few 
simple scenarios that make some definite predictions. In fact, as the number 
of technical assumptions in the models is reduced, we will find that the values 
of 612 and 613 point to a deviation from maximal 823 of the required size. 

The predictions of flavour models are meant to motivate the future exper- 
imental program, including precise neutrino oscillation experiments, as well 
as measurements of the absolute neutrino mass scale and of the neutrino-less 
double-beta (0u2j3) decay rates. In this spirit, we will consider lepton mass 
matrices that depend on a small number of parameters, and we will pro- 
pose a systematic procedure to illustrate the associated predictions (section 
2). We focus our attention on the mass matrix structures where the 2 — 3 
mixing angle is predicted close to the present best fit value, sin 2 623 — 0.4 
or sin 2 623 — 0.6. We will demonstrate that they can be associated to sim- 
ple, spontaneously broken, discrete flavour symmetries. One scenario cor- 
responds to neutrino mass matrices with two zero entries and two non-zero 
entries equal to each other (section 3). A second scenario resorts to the sum 
of two contributions to the neutrino mass matrix: a flavour off-diagonal 
term, plus a flavour-universal one (section 4). Finally, the minimal model 
with a flavour symmetry broken by a flavon doublet is presented (section 
5). These models allow, in particular, to correlate the 2 — 3 mixing octant 
with the value of the CP violating phase 5, as well as with the ordering of 
the neutrino mass spectrum (normal or inverse), thus motivating the exper- 
imentalists to solve the degeneracies that affect the long-baseline oscillation 
program [7]. 
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2 Prom the neutrino mass matrix to the 
observables 



We assume the Standard Model is extended to include a Majorana mass 
matrix for the three active neutrinos. The nine flavour observables in the 
neutrino sector can be currently divided in two groups. On the one hand, 
four physical parameters are precisely determined (with less than 10% un- 
certainty): 

p a = Am% 1 , lAmJj, 0i2, 0i3, (1) 

with the sign of Am|i still unknown. On the other hand, five parameters are 
at most weakly constrained. The mixing angle 023 can largely depart from 
the maximal value, 0.34 < sin 2 023 < 0.67 at 3a [2 [3]. The lightest neutrino 
mass miight can vary between zero and a few tenths of an electronvolt, the 
upper bound depending on the combination of cosmological data; we will 
assume a conservative limit, J2i m i ^ 0-5 eV at 95% C.L. [8j. The Dirac- 
type CP-violating phase 5 is unknown and can lie anywhere in the interval 
[0,27r), with some range of values already disfavoured at the la level by 
global fits of oscillation data [3]. Finally, the two Majorana-type CP- 
violating phases are unknown as well. The only experimentally accessible 
quantity sensitive to the latter two phases is the effective 0^2/3-decay mass 
parameter m ee , which is experimentally known to be smaller than a few 
tenths of an electronvolt; we will adopt the upper bound m ee < 0.38 eV at 
90% C.L., coming from the EXO-200 experiment [9]. The central objective 
of the present and future experimental program is the precise determination 
of these four observables (together with the sign of Am 2 ^): 

Xi = 023, "flight, S, m ee . (2) 

Flavour symmetry models may predict special structures of the neutrino 
mass matrix M v (in the basis where the charged lepton mass matrix M e 
is diagonal), that depend on few parameters, and therefore predict corre- 
lations between the observables. In particular, four real parameters of M v 
can be fixed to reproduce the observed values of the p a 's, up to their small 
experimental uncertainties. If M u depends on four (or less) parameters 
only, then one can derive a unique prediction for all the other observables, 
xf = xf(p a ), where the superscript denotes the sign of Am^. In this 
case the precise measurement of one x% can rule out the given matrix struc- 
ture. If M u depends on five parameters, then it is sufficient to specify the 
value of an additional observable Xj to predict the values of the other three, 
xf = xf(xj,p a ), i 7^ j. Such matrix structures can be tested in general by 
measuring with sufficient precision two Xj's. And so on and so forth. 

If the procedure to compute xf is carried out analytically, then it pro- 
vides the full set of predictions for a given form of M v . In general such 
procedure requires some tedious algebraic manipulations, but the advan- 
tage is to avoid a numerical sampling of the allowed range for the input 
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parameters, or an expansion in powers of some small parameter, with the 
associated ambiguities. 

Let us illustrate our method for a specific matrix structure that will be 
useful in the following, and that depends on five physical parameters: 

(a b c \ 
b d . (3) 
c d / 

The entries a, b, c and d are in general complex, but three phases can be 
absorbed in a redefinition of the three neutrino wavefunctions. Thus, there 
is a one-to-one correlation between the sought-after observables x^s. The six 
functions xf = xf(xj) for i ^ j are shown in Figures 1 and 2, for the case 
'+' (miight = rni, normal ordering) and '— ' (miight = "^3, inverse ordering), 
respectively. The procedure to derive xf (xj ) analytically is described in the 
Appendix, and it applies to any matrix M v with two zero entries. In principle 
one may find an analogue procedure for any structure of M v that depends 
on five parameters (see section [4] for a particularly interesting example) . 

As the matrix in Eq. ^ maintains the same structure under conjugation, 
for 5 f-> —5 the other observables are unchanged, therefore in Figures 1 and 
2 we need to display only the interval [0, tt] for 5. This matrix structure 
is also invariant under the exchange of the 2nd and 3rd rows and columns. 
This corresponds to the following exchange in the values of the observables: 

7T 

023 <-»■ 2 ~ ^23 and 5 <H> it - 5 . (4) 

Such relation can be clearly observed in the plots. In the case of inverse 
ordering, illustrated in Fig. the 2 — 3 mixing lies in the first (second) 
octant for cos 5 > (< 0). The curve 8(623, Pa) is well described, up to 
small corrections of order Am^/Aml, by the relation 

cos (5 tan 2^23 tan 2^12 sin ~ 1 . (5) 

The mass spectrum can be quasi-degenerate or hierarchical, with a lower 
bound 77J3 > 0.02 eV, and m ee ~ 771,3. The upper bound on ^i m « implies 
sin 2 (9 2 3 < 0.48 with cos 6 > 0.13, or sin 2 6» 23 > 0.52 with cos 5 < -0.13. In 
the case of normal ordering, illustrated in Fig.[2j the 2—3 mixing is extremely 
close to maximal, the mass spectrum is quasi-degenerate, with mi > 0.15 
eV. As a consequence, values of 5 close to or tt are predicted. The equality 
777 ee ~ 777j holds in very good approximation. The present upper bound 
on Y2i m i almost rules out the normal ordering case. Our exact results are 
consistent with the analysis of the matrix in Eq. ([3]) performed in Ref. [11] , 
where the correlations between the Xj's were understood by a sequence of 
careful approximations. 

In order to identify a possible flavour symmetry responsible for the mass 
matrix structure in Eq. (I3J), it is useful to inspect the size of the parameters 
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Figure 1: The correlations between sin 2 #23; S, ms (in eV) and m ee (in eV) 
for the neutrino mass matrix in Eq. Q), in the case of inverted ordering of 
the mass spectrum. The thick purple lines correspond to the best fit value of 
the parameters p a 's in Eq. while the purple shaded regions correspond 
to the 3a allowed range for the p a 's. The yellow bands are excluded by oscil- 
lation experiments. The green (blue) bands are excluded by the cosmological 
upper bound on ^ mi (by the EXO-200 upper bound on m ee ). These three 
exclusion bands are taken into account in all the six panels: the excluded 
portions of the best fit lines are dashed and the excluded portions of the 3a 
regions are not shaded. 
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a, b, c, d, that are needed to reproduce the experimental data. Note that, 
once the four input values for the p a 7 s are fixed, the mass matrix elements 
are fully determined as a function of one observable Xj. In Fig. 3 we show 
m a p = KMjy)^) as a function of sin 2 #23, for a(3 = ee,efi,er, [it. One 
observes that the pair m ee and m^ T is quasi-degenerate in the whole range 
of 023 , while the pair m eA1 and m eT becomes more and more degenerate as 
$23 approaches the maximal value. Note also that, in the case of inverse 
ordering, the two pairs cross each other for some special values of 623, that 
actually lie in the present la preferred regions. 
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Figure 3: The absolute values m a p (in eV) of the entries of the M u in 
Eq. (£5]), as a function of sin 2 623, for the case of inverse (left-hand panel) 
and normal (right-hand panel) ordering. The thin solid blue line corresponds 
to m ee , the thick dashed red line to mu T , the thin dashed blue line to m eil 
and the thin dotted blue line to m eT . The shaded bands are the same as in 
Fig. [7J The two pairs of vertical dashed lines limit the la preferred region 
for 623 in each octant. 



A symmetry for zero and equal entries in the 
neutrino mass matrix 



In the previous section we proposed a procedure to illustrate the phenomeno- 
logical implications of a neutrino mass matrix depending on few parameters. 
From a theoretical point of view, one shall wonder whether the minimal vi- 
able matrix structures can be justified by a symmetry. 

In this section we will focus on the matrix in Eq. ^ in the limit b = d, 
since it corresponds to two values of sin 2 623 close to the presently preferred 
values in the first and in the second octant, as shown in Fig. [3} Taking 



M„ 




(6) 



and considering the 3a allowed range for the parameters p 's, one finds two 
possible solutions: 



' sin 2 



23 = 0.40 
cos 5 = 0.59 



"Might = m 3 
m PP = 0.036 



+0.02 
0.01 
0.12 

0.14 

- m7+ - 001 pv 

- u.uof_ 002 ev 



' sin 2 623 
cos 5 = 
"Might 



+0.002 
^0.001 



eV 



0.62 



+0.03 
-0.02 



-0.75 



+0.15 
-0.12 



ni3 = 0.0289 



+0.0002 
-0.0001 



eV 



m,, 



0.02841°;°°°° eV 



(7) 



At first sight it may seem arduous to justify the vanishing of two entries 
of M v together with the equality of two other entries, at the same time 
maintaining the charged lepton mass matrix M e diagonal. Here we will 
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illustrate the ingredients needed to achieve this result and build a represen- 
tative flavour symmetry model. We remark that, no matter how elaborate 
such model may be, it provides the sharp predictions in Eq. ([7]) for the 
observables. 

The search for the appropriate flavour symmetry goes as follows: 

• The equality of the matrix entries M e „ and M T ^ requires that the two 
lepton doublets l e and l T transform together under the symmetry, that 
is, they belong to a dimension-two irreducible representation (irrep) 
2i, while Ifj, is a flavour singlet 1;. In this way the two equal matrix 
entries can arise from a flavon doublet (f), that is, a scalar field in a 
doublet irrep 2^, that acquires a vacuum expectation value (vev) in 
the direction (0) oc (1,1). 

• The tensor product 2; x 2; must lead to a non-zero diagonal entry, 
M ee = a, and a vanishing one, M TT = 0. Some simple group theory 
leads to the conclusion that this tensor product must contract with an 
additional flavon doublet cj)' , with (ft) oc (1,0). Since <p and <\>' must 
contribute separately in the sector 2\ x 1/ and 2/ x 2/, the flavour group 
must provide two different doublet irreps, 2^ ^ 2^. 

• Finally, since Man vanishes, the tensor product 1/ x 1; should not be 
invariant under the flavour symmetry. 

The simplest group that satisfies these properties is the order-twelve quater- 
nion group Qq, also known as the double-dihedral group D' 3 , that has two 
doublet irreps 2i and 22, plus four singlet irreps lj, with in particular 
I3 X I3 = I4 X I4 = I2, different from the invariant singlet lx- For the 
relevant group properties and conventions, we refer to Ref. [10], that pro- 
vides useful tools for the analysis of (double-)dihedral symmetries in general. 

Let us begin from the neutrino sector, writing the Majorana mass term 
as Uijlilj4>k(hh/ A), where h is the Higgs doublet and are the flavon fields 
that transform under Qq (we normalized (f>k to be dimensionless, that is, a 
scalar field over some cutoff scale). The lepton doublets and the flavons are 
assigned to the irreps of Qq in order to satisfy the above list of requirements, 



2 2 , l„~l 8 , (^)~2 2 . (8) 



JtJ \(pid 
Then, the Q6-invariant Lagrangian has the form 



y eT l e l T + ^y e e(lele<t>2u + Mr^d) + Ve^e$\u + Wld)^ 



hh 

— + h.c. 
A 



(9) 

Taking = ci(l,l) and (^2) = 02(1, 0), one obtains the neutrino mass 
matrix in Eq. with a = y ee C2V 2 /A, b = y efl civ 2 /A and c = y eT v 2 / A. 

The model works as long as M e is diagonal in the same basis. The 
charged lepton mass term can be written as X^Mej^kh* , where are the 
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charged lepton singlets and the flavons, in general different from those 
in the neutrino sector. With the assignment 

(::)— >* (:::)-„ ao, 

the (56-invariant Lagrangian reads 

- C e = [\ e (l e e c <pid ~ It^ C( P\u) + A At Z^(/uV2u + e c ip 2 d) 

+ X T (l T cp ld + l e ip lu )T c ]h* + h.c. . (11) 

For (ipi) = 0^(0,1) and (1^2) = c 2 (l,0), the charged lepton mass matrix is 
diagonal, with m e = X e d x v, = X^c' 2 v and m T = X T c^v. 

This completes the prove of existence of a minimal model leading to 
Eq. Note that the flavour group Qq is entirely broken by the flavon 
vevs, both in the neutrino and in the charged lepton sector. 



4 Off-diagonal plus universal neutrino masses 

Here we propose another, remarkably simple structure of the neutrino mass 
matrix that depends on five physical parameters and predicts the preferred, 
non-maximal #23. We will justify this structure with a minimal flavour 
symmetry at the end of this section. 

Our starting point is the so-called Zee mass matrix \12\ H~3] , with only 
off-diagonal entries, 

/0 a b\ 

Mf = la c , (12) 
\b c 0/ 

in the basis where M e is diagonal. The phases of a, b and c can be rotated 
away, thus there is no CP-violation and the leptonic mixing is described by 
a real orthogonal matrix R, with M° ff = R diag(Ai, A2, A3) R T . A simple 
way to extract a constraint on the mixing angles is to regard the three zero 
diagonal entries as a homogeneous system of three linear equations, 



0, i = 1,2,3, (13) 



where Zy = Rfj. In order to have a non-trivial solution, the matrix Z 
must be singular. The equation det Z = establishes a relation between the 
mixing angles, 



1 1 2(1 -2 sin 2 6 



13 J 



tan(20 12 ) tan(26» 23 ) sin6> 13 (l - 3 sin 



2 , 



'13 ) 



±1 , (14) 



where the positive (negative) sign corresponds to 623 < 7r/4 (623 > tt/4) 
and 5 = (5 = it). Not surprisingly, this equation has been already de- 
rived in various ways in the literature (see e.g. Refs. [HI [IS]). Indeed, the 
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Figure 4: The predicted value o/sin 2 #23 as a function o/sin 2 #12 and sin 2 #13, 
for the off-diagonal neutrino mass matrix of Eq. (12). The green shaded band 
corresponds to the la range for sin 2 6*23 in the first octant (left-hand panel) 
and in the second octant (right-hand panel). The vertical (horizontal) band 
corresponds to the 3a range for sin 2 #12 (sin 2 #13,). Dashed lines are the best 
fit values of the mixing angles. 



corresponding prediction for 622, fulfills with great accuracy the most recent 
oscillation data, as shown in Fig. [4} Taking the best fit values and the 3a 



ranges for 6*12 and #13 from Ref. [3], Eq. (14) implies 



sin 2 23 = O.42^ 2 or 0.58^2- (15) 

It is remarkable that these intervals roughly coincide with the la preferred 
regions for 623- 

The Zee mass matrix is known to be unable to reproduce all current data 
P~H EH H7j: it was shown that, if one fits the two neutrino mass squared 
differences, then the resulting mixing pattern is close to bi-maximal, which 
rules out the model. This is simply because the three parameters a, b and 
c are insufficient to accommodate the values of the four precisely measured 



observables p^s. In fact, the three zeroes on the diagonal imply Eq. (14) as 
well as two relations between the eigenvalues, 

Ai + A 2 + A 3 = 0, Ai(cos 2 #i2-tan 2 # 13 ) = A 2 (tan 2 13 - sin 2 12 ) . (16) 

Given the measured values of the mixing angles, the mass squared differences 
Am?. = A 2 — A 2 are incompatible with data. 

Now, the perfect prediction for 623 makes us wonder whether the lepton 
mixing pattern is indeed determined by an off-diagonal mass term, while 
the neutrino mass eigenvalues are not. Then, the most natural extension of 



Eq. (12) is the addition of a universal mass term, 

/0 a b\ (1 
M v =\a c + d 1 I , (17) 





a 






f 1 



















1 






c 


0) 




Vo 
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with d = \d\e lx complex in general. This modification does not affect the 
mixing matrix, as long as a, b and c are real (or have the same phase). The 
neutrino mass eigenvalues are given by 

mi = \\i + d\ , (18) 

and one is able to accommodate the two neutrino mass squared differences. 
In particular, we obtain a condition on the mixing angles that determines 
the ordering of the neutrino mass spectrum: 

normal for s\ 2 < 13 ^~ — , inverse for s\ 2 > 13 ^~ — • (19) 

The experimental values of 9\2 and #13 prefer normal ordering, with the 
inverse one only marginally allowed at 3<r. 



Note that the M u in Eq. (17) depends on five physical parameters, but it 
behaves in a singular manner with respect to our general analysis of section[2j 
the four input parameters p^s fix a, 6, c (or equivalently the three Aj's) and 
the product \d\ cosx- Two output parameters are automatically determined 



by Eq. (14): 623 and 5 (equal to or ir). As d is in general complex, there 
is CP-violation, but only in the form of Majorana phases. The other two 
output parameters, mught and m ee , depend on \d\ and cosx separately, and 
thus there is a non-trivial correlation between them, m ee = m ee (mn g h t ,p a ), 
that is illustrated in Fig. [5] Since the product \d\ cosx is fixed by the p a 's, 
there is a lower bound on m ee = \d\, which translates into a lower bound 
on mjight too. In the case of normal ordering, taking into account the 3a 
uncertainties on the PaS, we obtain 

mi > 4.3 x 1(T 2 eV, m ee > 2.1 x 1(T 2 eV. (20) 

The inverse ordering is disfavoured by the cosmological bound on ^ mi. 



Let us show that the matrix in Eq. (17) can be motivated by a minimal 
flavour symmetry. Clearly, the three lepton doublets shall transform in a 
triplet irrep of the flavour group. The smallest group with such an irrep is 
A4, the alternating group of 4 objects (we adopt the conventions in Ref. p~5| ) . 
The universal mass term is precisely the unique ^-invariant contribution 
to M u , while the off-diagonal mass term is generated by an ^-triplet flavon 
field (01,02,^3): 



+ y O ff(01^T + le&lr + lelfifo) 



hh , 
— + h.c. . 21 
A 



Note that, in order for a, b and c to have the same phase, the three vevs 
(4>i) should be real, that is the case as long as the flavon potential does 
not break CP. The other requirement is to keep M e diagonal, with three 
different eigenvalues. This is promptly obtained as soon as the charged 
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Figure 5: The effective mass m ee as a function of mught f or the neutrino 
mass matrix of Eq. (Tl\ ), in the case of normal (left-hand panel) and inverse 
(right-hand panel) ordering. The solid line corresponds to the best fit input 
parameters p a 's, while the purple shaded regions correspond to their 3a al- 
lowed range. The green (blue) band is disfavoured by the cosmological upper 
bound on ^rrij (by Ov2f3-decay searches). 



lepton isosinglets e\ also transform as an ^-triplet. Introducing flavons <p' 
and <p" in the irreps 1' and 1" of A±, respectively, one obtains 

-C e = \yi^2l a e c a + y v ^'{l e e c + uj1^ c + uj 2 I t t c ) 

a 

+ y vl <p"{l e e c + u\fj, c + loI t t c ) h* + h.c. , (22) 

where uj = exp(2-7ri/3). This leads to M e = diag(m e , m», m T ) as desired. 
Note that A4 is broken to nothing in the neutrino sector, because one needs 
( < H) (^2) 7^ (^3) 7^ 0, while in the charged lepton sector a Z2 X Z2 
subgroup is preserved. 



5 Minimal 1+2 lepton flavour symmetry 

In the sections [3] and [4] we started from predictive mass matrix structures 
and built a flavour model to justify them. Here we follow the opposite path: 
we will construct a flavour model, with the minimal possible number of 
technical ingredients, and study the corresponding correlations between the 
observables. 

Let us hypothesize that the three families of lepton doublets U as well as 
the three charged lepton singlets e| transform in a dim-one and a dim-two 
irreps of a flavour symmetry group: 

U ~ (1/ + 2,)< , A ~ (l e c + 2 e c)i . (23) 

In order to generate a non-zero mixing between the flavour singlet family 
and the doublet ones, it is mandatory to introduce a flavon in a dim-two 
irrep, <f> ~ 2$, that acquires a non-zero vev. In the search for a minimal 
theory of lepton flavour, it is meaningful to ask whether a viable model 
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exists, that contains only flavons of this type, and what are the associated 
predictions. 

The possible contributions to the lepton mass matrices are given by the 
irrep tensor products that are invariant under the flavour group: 



(lj x 


(1; x2jX 2^)i 




(2, x 2,)i + (2i x 2i x 2^)i 



m v ~ v ' t;i ,„ 7o \:;^o s , (24) 



(lj X leOl 


(1; X 2 e c X 20e)x 


(2/ x l e c x 2^)1 


(2, x 2 e c) a + (2, x 2 eC x 2 0e )i 



(25) 



where the subscript 1 denotes the invariant component in the tensor prod- 
ucts, if any. Note that we included at most one power of the flavon irrep; 
this is the whole story if the theory is renormalizable and the flavon fields 
carry also gauge indexes, that is to say, <\> v is a Higgs triplet, and cp e a 
Higgs doublet. Alternatively, one can treat the flavons as gauge singlets; in 
this case higher dimensional operators, suppressed by a cutoff scale A, may 
contain several powers of the flavons, leading to more complicated flavour 
structures; here we will neglect these corrections of order ({(j))/A) n . 

Let us identify the simplest viable mass matrices with the structure 



given in Eq. (24) and Eq. (25). The minimal possibility is to employ a 
unique dim- two irrep (2j = 2 e c = 2$ = 2), and to assume that the tensor 
product 2x2 contains a singlet, while 2x2x2 does not (to avoid an 
additional, unnecessary invariant). This is the case for the two non-abelian 
order-eight groups -D4 and Q = D' 2 , that satisfy 2x2 = Ylt=i ^ii where li is 
the invariant singlet (we use conventions as in Ref. |18| ) . In the case of the 
quaternion group Q, the invariant component in 2\ x 2; is antisymmetric, 
therefore the corresponding contribution to M v vanishes; we will see that 
such contribution is needed to accommodate the data, thus we will rather 
focus on the dihedral symmetry D4. In this case (2; x 2/)i = + fah and 
the lepton mass matrices take the form 

a b 

M u =\ ... c 1, M e = I D C I , (26) 




where the diagonal entries are independent from the flavons, while the off- 
diagonal ones are linear in the flavons. According to our minimality tenet, we 
put to zero the 11-entry of M e by choosing 1; x l e c ^ 1^ Denoting the flavon 
vevs by (P) = (P u r d ) T and <c/> e ) = (^ ^) T , one has b/b' = ±(<%/<%) ± \ 
B/B' = ±(^/0|) ±1 and D/D' = ±(^/^) ±1 , where the signs are not 
correlated and depend on the choice of the irreps 1; and l e c 
The physical mixing matrix Upmns is determined as usual, 

M v = U* u d v Ul , M e M] = U* e d 2 e Uj , Upmns = U\U V , (27) 

where d u and d e are diagonal matrices containing the neutrino and charged 
lepton mass eigenvalues. Let us remark that the vev of a doublet 4> breaks 



13 



-D4 completely, unless it is aligned in the directions (1,±1), (1,0) or (0, 1): 
in these cases a Z<i subgroup is preserved. We will consider only these vev 
alignments, that are protected by the residual symmetry and correspond to 
a minimal number of parameters in the mass matrices. In the case 6 = 
(or V = 0), My is diagonalized by a single rotation of angle 9 U ; the case 
b = ±6' can be reduced to the previous one by a it/ 4 rotation in the flavour 
doublet sector. Similarly, for B = D = (or B = D' = 0, etc.), M e M} 
diagonalized by a unique rotation 9 e , while for B = ±B' and D 



is 

±D' an 

additional maximal rotation is needed. If this maximal mixing is present 
(or absent) in both M v and M e Mj, it cancels out in Upmns and it is 
straightforward to check that one entry of Upmns vanishes, at odds with 
experimental data [3j. Thus, data indicate that <f> v and <f> e should be two 
different flavons, one acquiring a vev in the direction (1, 0) or (0, 1), and the 
other in the direction (1,±1). 

The general idea to break a flavour symmetry group G to two different 
subgroups G u and G e in the neutrino and in the charged lepton sector, re- 
spectively, is not new: it was mostly employed to explain the tri-bi-maximal 
mixing pattern, in models where the three lepton doublets transform in a 
triplet representation of G [19]. Here we have shown that, in the context of 
1 + 2 flavour models, the minimal viable possibility is to take G u = Z2 and 
G e = Z 2 . 

Let us derive the predictions of this scenario. Following the above con- 
siderations, in Eq. (26) we take b' = 0, B' = B and D' = —D (there are a 
few other equivalent assignments). Then, one can check that M v and M e M\ 
depend on 6 independent absolute values of the matrix entries, plus 2 phys- 
ical phases, therefore the 6 lepton masses, the 3 mixing angles and the 3 
CP-violating phases are correlated, as we now describe. One finds 



cos 9 u e %a sin 9 u e 
U v = I — sin 9 u e ia cos 9 v e 
1 




d v = diag(m_, c, m+) , (28) 



where a = arg(a6* + bc*)/2, tan 2^ = 2\ab* + bc*\/(\c\ 2 - |a| 2 ), |m ± | 2 = 
\b\ 2 + \a\ 2 /2 + \c\ 2 /2 ± ^{\c\ 2 - \a\ 2 ) 2 /4 + \ab* + bc*\ 2 . Note that the two 
independent mass squared differences Am^ can be accommodated indepen- 
dently from the value of two parameters a and 9 V , that enter in Upmns- in 
the charged lepton sector one has 



U P 



\ 



cos 9 e 


sin 9 e e 1/3 





\ 


sin 9 e e l13 


cos 8 e 


1 




V2. R 


V2 


~V2 




sin 9 e e l P 


cos 9 e 


1 




V2 


V2 


71 


/ 



d e = diag(0,TO^,m r ) , (29) 



with (3 = arg(5C*), m 2 = 2\B\ 2 + \C\ 2 , m 2 = 2\D\ 2 + \C\ 2 and tan# e = 
\f2\BjC\. Note that m e vanishes, which is a very good approximation, and 
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Figure 6: The correlation between sin 2 823 and 5 defined by Eq. (30). The 
thick purple line (the shaded purple region) corresponds to the best fit values 
(the 3a allowed ranges) for 812 and 813. The vertical yellow bands (dashed 
lines) limit the 3a (la) allowed range for sin 2 823. 



the two masses and m T are independent from the two parameters 8 e and 
(3, that enter in Upmns (if the assignment of and m T were interchanged, 
6 e would be forced to lie close to ir/2 and the observed mixing angles could 
not be reproduced). 



By using Eq. (27), one can easily check that 612, #13, #23 and 8 are a 
function of 8 U , 9 e and a unique combination of phases, 2a — /?, while the 
two physical Majorana-type phases are a function of the other independent 
combination of phases too. The correlation between the mixing angles and 
5 depends on the assignment of the neutrino masses: we discuss first the 
case mi = l^-l) m 2 = \c\ and 7713 = \m + \, that corresponds to normal 
mass ordering. Comparing the predicted form of Upmns with its usual 
parameterization, we obtain 

2s 12 c 23 s 13 + 4S12C12S23C23S13 COS 5 + 2c 2 2 S23 = 1 • 

This relation is illustrated in Fig. [6j For the best fit values for 812 and 
013, one finds cos 5 > 0.49, and for 5 = one reaches the minimal value 
sin 2 #23 — 0.616, on the edge of the la preferred region in the second octant. 

Adding also Am^ and Am^ as input parameters, one can determine 
the lightest neutrino mass mi and the 0^2/3-decay effective mass m ee as 
a function of 5 (or equivalently #23) an d of the extra phase a. This is 
because, in this scenario, in the basis where M e is diagonal M u depends on 
six parameters, therefore, after the four jj a 's are fixed, one can still vary two 
independent CP-violating phases. One can derive the explicit dependence of 
mi and m ee on a, by generalizing our procedure of section [2] For simplicity 
we limit ourselves to the CP-conserving case, because Eq. po] ) favours a 
value of 5 close to zero, and because the value of m ee for a CP-violating 
value of a is intermediate between its two values in the CP-conserving case: 
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for 5 = we obtain mi ~ 0.036 eV and, for a = (a = vr/2), m ee ~ 0.012 
eV (m ee ~ 0.034 eV). 

We remark that Eq. (30) follows from our minimality assumptions, that 
is to say, it is the consequence of the simplest realization of a 1+2 flavour 
symmetry. This implies a non-maximal value of 623 in the second octant. 
It is interesting to compare Eq. (30) with the relations among the mixing 
angles Oij's and the phase 5, that have been obtained in Ref. [6]. In that 
approach one assumes that the flavour symmetry group is generated by a Z2 
symmetry in the neutrino sector plus a Z m symmetry in the charged lepton 
sector, with m > 3, and this leads to two constraints on the values of 9%j's 
and 5. The same analysis can be carried on in the present case, with m = 2, 
and we find that one of the two constraints is trivially satisfied, while the 
other is given by Eq. (30). 

If the assignment of the neutrino mass eigenvalues is mi = |m_|, rri2 = 
\m+\ and m.3 = |c|, both normal and inverse orderings are possible, and 
Eq. ( 30 ) is replaced by 

= 1 . (31) 

In this case the best fit value of #13 corresponds to an almost maximal 2 — 3 
mixing, with sin 2 623 — 0.49, a value (slightly) disfavoured by present data. 
The other assignments of the three neutrino mass eigenvalues are not viable. 



6 Conclusions 

We attempted to provide a modern view on the lepton mixing, and to pro- 
pose some keys to interpret it. The underlying flavour dynamics may be 
intricate and difficult to disentangle, especially if it occurs at very large en- 
ergy scales. Nonetheless, the effective structure of the lepton mass matrices 
may well be simple and depend on a small number of parameters. As a 
matter of fact, this must be the case for a flavour model to be testable. 

Thus, it is useful to identify the minimal lepton mass matrices M v and 
M e compatible with those observables that have been measured precisely: 
the charged lepton masses, the mixing angles #12 and #13, and the neutrino 
mass squared differences Aro| . Such viable flavour structures imply definite 
predictions for the neutrino mass observables that are not precisely measured 
yet: the mixing angle 623, the CP-violating phase S, the lightest neutrino 
mass miight and the 0^2/3-decay effective mass m ee . We proposed a general 
procedure to illustrate the correlations between these observables. 

In particular, it is intriguing that several minimal flavour structures im- 
ply a deviation from maximal 2 — 3 mixing of the size that is presently 
suggested by the data. We identified a few classes of mass matrices with 
this property: 

• A matrix M u with two zero entries and two other entries equal to each 
other, in the basis where M e is diagonal. There are two remarkable 
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Figure 7: The interval for sin 2 623 displayed here is the 3a allowed range 
from a global fit of oscillation experiments J3j/. The lowest, blue dots (bands) 
are the central values (the la ranges) in the two octants. The upper, red 
dots (bands) are the predictions of our models, assuming the central values 
(the 3a ranges) for the four parameters in Eq. |7]). The three predictions 
correspond to Eq. ([?]] (Q§ model), Eq. (A4 model), and Eq. (30) with 

5 = (D4 model with no CP violation). 



cases: [il M. 



and M, 



e/i(er) 



M Mr or M, 



efj,(er) 



M P , 



(section [2]); [ii] M ee = M e ^ (eT) = and M^ T = M jU(U ( Tr ) (appendix). 
We pointed out the properties of the flavour symmetry group that is 
needed to explain these matrix structures (section [3]). 

• An off-diagonal contribution to M u plus a contribution proportional 
to the identity, in the basis where M e is diagonal. Such structure is 
tailor-made for an A4 flavour symmetry (section [4|. 

• Matrices M u and M e that respect each a different Z2 symmetry, in 
such a way that only one free mixing angle is present in each sector. 
The minimal realization is provided by a D4 symmetry, with the three 
lepton families transforming in the 1 + 2 irreducible representations 
(section [5]). 

The three specific models that we presented predict 623 close to its best fit 
values, as summarized in Fig.[7j Thus, the future precision measurements of 
this mixing angle will be able to confirm or rule out these scenarios. In order 
to discriminate among them, one shall use the correlated predictions for the 
mass ordering as well as for 5, mught and m ee , that could be all precisely 
measured within the next decade. 
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Appendix 

Here we detail the procedure to derive the functions Xi(xj,p a ), defined in 
section [2j for a special class of Major ana neutrino mass matrices M v that 
depend on five physical parameters: the 3x3 symmetric matrices with two 
independent entries equal to zero. There are 15 such matrix structures, 
and only 7 of them are compatible with oscillation data, as first shown in 
Ref. [20] (see Ref. [21] for a recent analysis). 

The four constraints relating the observables (3 neutrino mass eigenval- 
ues, 3 mixing angles and 3 phases) have the form 

= M a /3 = y^ u {U*PMNs)ai{U*PMNs)l3i m i , (32) 

i 

for two independent complex entries M ctl p 1 and M a2 p 2 . Since these equa- 
tions are linear in rrii, they can be rewritten in the form 

— e^ 13 = fWij,S) , — = r 23 (9ij,S) , (33) 
m 3 1713 

where (fij are the Majorana phases and the functions depend on the 
choice of {ax, fix) and {a 2 , (3 2 ), 

^ _ Uai3Up 1 3Uct2lUp2l — Ug 1 lUl3 1 lUa 2 3Ul323 
UailUfaxUatfUptf — UatfUptfUaziUfax 

where U a i are the Upmns entries without the Majorana phases. The neu- 
trino mass squared differences can be expressed in terms of the neutrino 
masses as 

f Am 2 21 =m 2 {\f 23 \ 2 -\f 13 > 2 ' 
\ Am 2 31 = m|(l - |/ 13 

These equations are linear in m 3 and linear or quadratic in cos 5. Therefore 
the system has at most 2 solutions (for each given ordering of the mass 
spectrum) of the form 

ml = ml(6 23 ,p a ), cos 5 = cos5{9 23 ,p a ), (37) 



2 .2,, |, ,2, ( 36 ) 
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Figure 8: The CP-violating Dirac phase 5 (top left), the lightest neutrino 
mass mi in eV (top right) and the absolute values of the neutrino mass 
matrix elements m a p in eV (bottom), as a function of sin 2 #23, in the case 
m ee = m efl = 0. The ordering of the mass spectrum is necessarily normal. 
The yellow bands are excluded by the 3<r bound on 623; the vertical dashed 
lines limit the la preferred region for 623- In the bottom panel, the thin 
dotted blue line corresponds to m er , the thick dotted red line to m^ T , the 
thick dashed red line to m aa and the thick solid red line to m TT . 



where the input parameters p a 's are defined in Eq. ([T]). Each solution is 
physical only when it fulfills the requirements m\ > and | cos<5| < 1. The 
effective Oz/2/3-decay mass parameter can be written as 

mee = m 3 I E& 2 /*3 + U*Jf23 + Kl \ = m ee (023, P a) , (38) 

where one has inserted the expressions for 777,3 and S given by Eq. p7| ). 
Analogously, any matrix element M a p can be calculated as a function of 623 
and the P aS. The functions Xi(023, Pa ) for the observables xi = 777,3,5, m ee 
provide a parametric representation for the functions Xi(xj,p a ), for any pair 
This concludes the analytic computation of the observables for any M v 
with two zero entries. 

The results for M M(U = M TT = (case C) were presented in section [2j 
Here we present the results for the 6 other viable cases: 



Ai 




= M ee 


= 0, 


A 2 : 


M er 


= M ee 


= 0, 


(39) 


Bi : 


M eT 


= 


= 0, 


B 2 : 




= M TT 


= 0, 


(40) 


B3 : 




= Mw 


= 0, 


B A : 


M eT 


= M TT 


= 0. 


(41) 
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The matrices in the second column are obtained from those in the first col- 
umn, by the exchange of the 2nd and 3rd rows and columns, therefore their 
predictions are simply related by Eq. Q. Thus, we will illustrate the results 
for the cases A\, B\ and B3 only, that for normal mass ordering correspond 
to 623 in the first octant, as preferred by global fits [21 13]. When inferring 
the results for the 'specular' cases, one should recall that the allowed exper- 
imental range for 623 is slightly asymmetric between the first and the second 
octant. 

For simplicity we fix the input parameters p a 's to their best fit values: the 
predictions are slightly different when the small uncertainties on these pa- 
rameters are taken into account. We present the plots for mii g ht(sin 2 623, p a ), 
5 (sin 2 623, p a ) and m a p(sm 2 623, p a ), that are sufficient to infer the correla- 
tion between any pair of observables x^s. 

• The results for the case A± are illustrated in Fig. [8j 

— Only a normal ordering of the neutrino mass spectrum is allowed. 

— The 2 — 3 mixing satisfies sin 2 823 < 0.505, however this upper 
bounds relaxes substantially when allowing the p a 's to vary in 
their 3a allowed range. 

— The phase 5 satisfies cos 5 < 0.62. 

— The lightest mass mi lies in the range [3.1, 6.3] x 10 -3 eV. 

— There is no contribution to the 0^2/3-decay rate. 

— In the la preferred region for 623, the equality m w = m^ T can 
be realized. Therefore, as in the case C discussed in section [2j 
the vanishing of two matrix entries and the equality of two others 
imply the desired prediction for 623. One could build a model on 
the lines of the one in section [3j in order to justify this matrix 
structure by a discrete family symmetry. 

• The results for the case B\ are illustrated in Fig. [9j 

— Both normal and inverse ordering are allowed. 

— The angle 623 lies in the first (second) in the case of normal 
(inverse) ordering. This statement remains valid when allowing 
the paS to vary in their 3a allowed range. 

— The phase 5 is very close to tt/2 (maximal CP- violation), with 
I cos<5| < 0.13 (a similar bound holds when the 3a uncertainties 
on the p a 's are taken into account). 

— There is a lower bound on the lightest neutrino mass, mi > 2.8 x 
10~ 2 eV for normal and 7713 > 2.6 x 10 -2 eV for inverse ordering. 

— There is a lower bound on the 0^2/3-decay effective mass, m ee > 
2.8 x 10~ 2 eV for normal or m ee > 5.5 x 10 -2 eV for inverse 
ordering. 
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Figure 9: The same as in Fig. [#| but for the case = m er = 0. The 
left-hand (right-hand) plots correspond to normal (inverse) ordering. The 
green (blue) bands are excluded by the cosmological upper bound on 
(the experimental upper bound on m ee ). In the bottom panels, the thin solid 
blue line corresponds to m ee , the thin dashed blue line to m eiI , the thick 
dotted red line to mn T and the thick solid red line to m TT . 



• The results for the case B3 are illustrated in Fig. 10. They are qual- 
itatively the same as in the case B\, with minor modifications of the 
predicted range for the Xj's. 



In summary, our procedure allows to illustrate clearly all the predictions 
that follow from the assumption of two zero entries in M v . We confirm that 
8 matrix structures of this type are already excluded by present data. The 
7 allowed matrix structures, analyzed above, could be ruled out by future 
measurements of the observables Xj's. No matrix with three zero entries 
is allowed, because there is no intersection between the predictions of the 
different matrices with two zero entries. 
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Figure 10: The same as in Figs. [#| and[$| but for the case m MA1 = m e ^ = 0. 
In the bottom panels, the thin solid blue line corresponds to m ee , the thin 
dotted blue line to m eT , the thick dotted red line to m^ T and the thick solid 
red line to m TT . 
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